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We have examined the effects of a bound dineutron, 2 n, upon big bang nucleosynthesis (BBN) as 
a function of its binding energy B2 n . We find a weakly bound dineutron has little impact but as £?2 n 
increases its presence begins to alter the flow of free nucleons to helium-4. Due to this disruption, 
and in the absence of changes to other binding energies or fundamental constants, BBN sets a 
reliable upper limit of P>2 n < 2.5 MeV in order to maintain the agreement with the observations of 
the primordial helium-4 mass fraction and D/H abundance. 
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PACS numbers: 99.99 



I. INTRODUCTION 



The concordance of the predicted synthesis of the 
lightest nuclei during the period immediately following 
the Big Bang with the observed primordial abundances 
presents us with a powerful probe of the state of the Uni- 
verse during its earliest epochs. In addition to constrain- 
ing standard cosmological parameters such as the density 
of baryons, the number of neutrino flavors and their de- 
generacy H S g 1 ft g 1 M El El El III E3 
BBN is sufficiently complex that one can also learn 
about such possibilities as Quintessence pH HH liH, fl9| , 
modifications of gravity fl3. Infill o r neutrino oscilla- 
tions /mass/decay Hi llfSllIl |H |H |H ^ • 

Perhaps the most intriguing use of BBN is in con- 
strainin g t he variation of the fundamental constants 
[H IH M HI IH HI . Support for this hypothesis 
has emerged from recent observations of quasar absorp- 
tion lines at redshift of z=l-2 by Webb et al. jH |3|| 
that suggest the fine structure constant, a, may have 
been smaller in the past (though see 39] ). In some cases 
the variation of a fundamental constant is easily imple- 
mented in BBN because the nuclear physics aspects of 
the calculation are unaffected, but in others the lack of 
an adequate theory to predict such parameters as the 
nuclear binding energies and cross sections introduces 
a degree of uncertainty. An example is the calculation 
of Kneller & McLaughlin [36j who looked at the vari- 
ation of Aqcd, its effects upon the binding energy of 
deuterons and the neutron-proton mass difference and, 
in turn, the impact upon BBN. In addition to variation 
of these nuclear parameters, variation of the constants 
relevant to nuclear structure might also partially stabi- 
lize nuclei that are presently particle unstable. For ex- 
ample, the lack of stable A = 5 and A = 8 nuclei is 
often cited as the explanation for the dearth of nuclei 
formed with masses above helium-4 though the endother- 
micity of pure strong reactions such as 4 He(T, n) 6 Li, 
4 He( 4 He, p) 7 Li and 4 He( 4 He, n) 7 Be plays a role. 



The focus in this paper is upon another nucleus that 
could also become stabilized - the dineutron 2 n. The 
dineutron, a member of the nucleon-nucleon isospin 
triplet, is a spin singlet and, by itself, the dineutron is 
weakly unbound 1 by ~ 70 keV, the n-n scattering length 
being negative [H EHH HI • Early direct searches [5| 
did not see evidence for a stable dineutron but recently 
Bochkarev et al. claim that 45 ± 10% of the decay 
of an excited state of 6 He is through the dineutron state 
2 and Seth & Parker [43, amongst others, find evidence 
for dineutrons in 5 H, 6 H and 8 He decay. There is also a 
claim for tetraneutron, 4 n, emission in the decay of 14 Be 

Given that the dineutron is only weakly unbound, even 
small changes in the pion mass could, perhaps, result 
in a bound dineutron, although at present there is not 
enough experimental information to show whether or not 
this would occur m |5(j ■ The scattering length is quite 
sensitive to the pion mass and so it is small changes in 
fundamental constants that change its mass, such as a, 
^QCD or the Higgs vev, that could cause the dineutron to 
become bound. Since we are lacking an exact relationship 
between these fundamental constants and the binding en- 
ergy of the dineutron we do not adopt a particular model 
for the time variation of fundamental constants, but in- 
stead explore the effect upon BBN of a bound dineutron 
directly. 

In this paper we shall make an effort to derive a con- 
straint upon the dineutron binding energy. We will con- 
sider the dineutron in isolation i.e. whatever the source 
of the new stability of the dineutron we shall limit the 
effect to just this nucleus. We begin with an overview 
of standard BBN in section ^H]with an emphasis on the 
details of the flow from free nucleons to hclium-4 before 
proceeding to insert dineutrons in section mill In section 
£|IVI we present our results for a baryon-to-photon ratio 
of 77 = 6.14 x 10~ 10 and follow it up in Jv]with a discus- 
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1 Though it may become stable on the surface of neutron-rich nu- 
clei. 

2 Bochkarev et al. |4(I also have evidence of diproton emission 
from an excited state of 6 Be. 
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sion of the errors in the calculation and any degeneracy 
with r\ in VII Finally, in HV11I wc show how BBN can 
limit the dineutron binding energy before presenting our 
conclusions. 



II. BBN WITHOUT DINEUTRONS. 

BBN can be simplistically broken into three phases 
characterized by the degree of equilibrium within the nu- 
cleons/nuceli. 

During the first phase, at temperatures above T > 
1 MeV, there are virtually no complex nuclei so that all 
the nucleons exist in a free state. The rapidity of the 
weak interactions in converting neutrons and protons, 



n «-> p + e + v e 
n + e <-> p + v e , 
n + v e «-> p + e, 



(la) 
(lb) 
(1c) 



establish a weak equilibrium so that the neutron/proton 
ratio, F, is simply F w exp(— A np /T). As the Universe 
cools the rate at which neutrons must be converted to 
protons in order to maintain the equilibrium cannot be 
accommodated. As a consequence, the neutron-to-proton 
ratio departs from its equilibrium value and is said to 
'freeze-out' even though conversion continues to occur. In 
the absence of neutron decay and the formation of com- 
plex nuclei, the ratio would attain an asymptotic value 
of F ~ 1/6 [66l |6tJ • The departure of F from its equilib- 
rium value denotes the boundary between the first two 
phases of BBN. 

During the second phase of BBN, from a temperature 
of ~ 1 MeV to ~ 100 keV, the abundances of the var- 
ious nuclei also begin to depart from equilibrium. At 
~ 1 MeV their abundances are suppressed relative to 
the free nucleons but the nuclear reactions that form 
them establish, and maintain, chemical/nuclear statisti- 
cal equilibrium (NSE). In equilibrium the abundance 3 , 
Ya = tia/ub, of a complex nuclei A is derived from 
fiA = Z fjrp + (A — Z) fi n so after inserting the expres- 
sions for the Boltmann number density we find 



for the thermal factors in equation J5J) we obtain 

A 3/2 



Y A 



9A^ ■ yl+Z-A yl-Z yA-1 

2[3V2F- 1 P n ° 



x exp 



B A -(A- 1)B L 
T 



(3) 



This equation now makes it much clearer that the abun- 
dance of a nucleus with mass A + 1 is suppressed by 
approximately Yd relative to the abundance of a nucleus 
with mass A. 

If the neutron/proton abundances are held fixed then 
the NSE abundance has a minimum at Ta = 2 Ba/(3A — 
3). Below Ta the various nuclear reactions provide suf- 
ficient nuclei to keep Ya in equilibrium as the abun- 
dance climbs from the minimum but eventually a point 
is reached where this required production rate cannot 
be met and the abundance falls beneath the equilibrium 
value. This departure from equilibrium for the complex 
nuclei is in contrast with that of the neutrons where it 
was an insufficient rate that led to the departure, here it 
is a lack of reactants that is the cause. Heavier nuclei are 
the first to depart from equilibrium: helium-4 departs 
at T ^ 600 keV while helium-3 and tritium drop out at 
T ~ 200 keV. Below T ~ 200 keV the only compound 
nucleus in NSE is the deuteron and its abundance con- 
trols the rate at which all the heavier nuclei can be pro- 
duced. By T ~ 100 keV the D abundance is approaching 
that of the free nucleons and the amount of D destruc- 
tion, via such reactions as D(D,p)T and D(D, n) 3 He, 
has becomes significant. When this occurs the deuteron 
abundance cannot be replenished sufficiently quickly to 
maintain its equilibrium and, consequently, it too departs 
from NSE. This final NSE departure forms the entrance 
to the third stage of BBN proper. 

Below the deuteron NSE departure temperature the 
D abundance continues to grow for a short period but 
eventually the D + D drain tips the balance in favor of 
destruction and the deuteron abundance reaches a peak 
amplitude. The tritons and helions formed via the reac- 
tions D(D,p)T and D(D, n) 3 He are produced in roughly 
equal amounts but the helions rapidly transform to tri- 
tium via 3 He(n,p)T. The last step in the formation of 
the alpha particle is almost exclusively T(D, n) 4 He which 
destroys ~ 1/3 of all the deuterons formed. The essential 
steps in the scheme are illustrated in figure Q). 
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Using F ~ 1/6 and a baryon- photon ratio of 77 ~ 10 10 
we see that for a temperature of T ~ 1 MeV the abun- 
dance of deuterons is Yd ~ 10~ 12 . After substituting Yd 



The term 'abundance' is also used for the ratio Ya/Yh 



Smith, Kawano and Malaney |&JJ identified 8 reactions 
among the A < 4 nuclei as being important for BBN. Five 
are identical to those in figure |JT|I. the sixth is neutron- 
proton interconversion, and the two reactions they in- 
cluded, and which we have omitted from the figure, are 
3 He(D, p) 4 He and T(p, 7 ) 4 He both of which are 2 - 4 or- 
ders of magnitude smaller than T(D, n) 4 He, so there im- 
portance is marginal. Whatever the exact number, this 
handful (or two) of important reactions is much smaller 
than the number of reactions included in any BBN code. 



3 



n 


+ 


P 










D 


+ 


D 



D 



+ 



T 



4 He 



3 He 



FIG. 1: A diagrammatic flow of nuclei in standard BBN. The 
complex nuclei are immersed in a (rapidly drained) bath of 
free nucleons and so we outline only where they form part of 
the set of reactants. 



tain an equilibrium and, typically, the flow of nuclei in 
any given reaction is one direction only. In a few cases 
where the Q- value for the reaction is less than few times 
the temperature during BBN (i.e. T ~ 100 keV, Q ~ 
500 keV) flow may be in both directions because the 
'activation energy' for an endothermic reaction is readily 
available. Examples of this behavior were seen in Kneller 
and McLaughlin |36( • The reaction should also be prefer- 
ably strong in nature since this is the behavior seen in 
standard BBN where reactions such as D(D, n) 3 He dom- 
inate over D(p, 7) 3 He though, in a few cases, such as 
p(n, 7)D, electromagnetic or weak interactions play an 
important role. The last requirement will remove those 
cases that pass the first two but whose entrance channels 
involve multiple particles. 

From these requirements we have selected four dineu- 
tron reactions that we expect to be important: 

p + 2 n^D + n, Q = 2.22 MeV - Bi n , (4) 
D + 2 n^T + n, Q = 6.26 MeV - B* u , (5) 



As the Universe cools eventually the Coulomb barriers 
in the various reactions become insurmountable leading 
to a cessation of the nucleosynthesis. The abundances 
have plateaued to their 'primordial' values with virtu- 
ally every neutron now incorporated in helium-4 and only 
small residues in D, T and 3 He. A small (but detectable) 
abundance of lithium-7 and beryllium-7 have also been 
formed but we will not discuss these two nuclei further. 



III. INSERTING DINEUTRONS INTO BBN 

Inclusion of a new, light nucleus into BBN has the po- 
tential to significantly influence BBN and alter the pre- 
dicted primordial abundances. These changes will oc- 
cur because the dineutron will disrupt the flow of nu- 
cleons through the reaction network by both presenting 
new exit channels for reactions that already exist in BBN 
and through new entrance channels in the formation of 
the nuclei. One can construct a large number of plausi- 
ble reactions in which dineutrons participate but not all 
are expected to play a prominent role for the same reason 
that standard BBN is dominated by only a few reactions. 
We can use our understanding of the important reactions 
in standard BBN to pick from the plethora of possibil- 
ities for dineutron reactions those which we expect to 
be important. The most important reactions involving 
dineutrons should be 

• preferably exothermic, 

• dominated by the strong interaction, 

• and two-bodied in their entrance channel. 



3 Hc + 2 n ^ 4 He + n, Q = 28.29 MeV - B* n , (6a) 
3 Hc + 2 n^T + D, Q = 2.99 MeV - Ba n . (6b) 

The dineutron binding energy will determine the Q- value 
in these reactions and, if we permit values of i?2 n of sev- 
eral MeV, both Q and (j6b|) may reverse sign. We will 
not consider any changes to the other nuclear binding 
energies. Though one may expect a large change in the 
dineutron binding energy to be reflected in equally signifi- 
cant changes to the structure of the deuteron, the effects 
upon three nucleon nuclei may be considerably smaller 
[52j . In our study, we also do not consider reactions lead- 
ing to the formation of destruction of nuclei above mass 
4. 

To this list we add three additional reactions: 

n + n^ 2 n + 7, Q = B2 n , (7) 
2 n <-> D, (8) 
2 n + p^T + 7 , Q = 8.48 MeV — £?2 n . (9) 

These reactions could become important for producing 
and, more importantly, removing dineutrons. In partic- 
ular, the inclusion of © and © is based on the follow- 
ing reasoning. Once we have dineutron cross sections 
we must integrate them over a Maxwell-Boltzmann spec- 
trum to obtain a thermally averaged rate [53 ■ In all the 
dineutron reactions the lack of a Coulomb barrier means 
that the cross section for an exothermic reaction varies as 
1/VE at low energy and so the rate becomes a constant. 
This can have important consequences: written in terms 
of the temperature, a reaction such as A + B — > X, with 
a rate per particle pair T, destroys nucleus A at a, rate 



Exothermicity plays a pivotal role in the BBN because, 
during its second two phases, the system does not at- 



(10) 
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where we have used the relations jib oc T 3 and T 2 oc 1/t. 
If we assume the abundance of nucleus B is much larger 
than A's and does not change by any other process 
then the solution to this equation is Ya oc exp(rYeo^ 1 ) 
where Ybo is the abundance of B at some fiducial tem- 
perature. In this scenario, the abundance of A never 
becomes a constant and BBN would never end! In 
standard BBN this situation never arises because the 
two temperature-independent reactions, p(n,7)D and 
3 He(n,p)T, are killed by the decay of the neutron. But 
if the dincutron becomes stable then without a reaction 
such as 2 n(p, y)T the dineutron abundance could plateau 
to a constant larger than, say, the abundance of helium-3 
and the circumstances of equation Hl()|) may be realized in 
the reaction 3 He( 2 n, n) 4 He. In order to obtain a primor- 
dial abundance of 3 He it is imperative that this situation 
be prevented. The decay of the dineutron and 2 n(p, 7)T 
will ensure this by depleting the final dineutron abun- 
dance for all values of i?2 n we shall explore. 

Now that we have determined the most important re- 
actions, we need cross sections for them before we can 
proceed. 



2 n(p,n)D, D( 2 n, n)T, 3 He( 2 n, n) 4 He and 3 He( 2 n,D)T 

Since the dineutron is presently unstable we posit cross 
sections based on 'similar' strong reactions involving deu- 
terium or other light nuclei. If we consider an arbitrary 
two-body strong reaction i+j <-* k+l then general consid- 
erations lead us to expect a cross section per particle pair 
which is proportional to a matrix element squared, the 
phase space with an energy conserving delta function and 
inversely proportional to a flux. In nuclear astrophysics, 
one typically writes the cross section as 



S(E) 
E 



exp 



-7T OL Z; Zr. 



2 ti-ij 
E 



(11) 



where /itjj is the reduced mass for incoming particles i 
and j, and S(E) is the astrophysical S-factor. For our 
purposes, this parameterization is not sufficient since it 
does not explicitly show the effects of a Q value upon the 
final states. This is particularly crucial since in our study 
Q-values will vary as the dineutron binding energy varies. 
With that in mind, we write the non-resonant (S-wave) 
contributions to the cross section, following EH , as 
proportional to the product of both Coulomb penetrabil- 
ity factors Gij(E) and Gki(E + Q), the available phase 
space in the exit channel $ki(E + Q) together with the 
statistical weight gu and the reciprocal of the entrance 
channel velocity. The penetrability factor for charged 
particle interactions, Gij(E), is simply 




exp 



— 7r a Z;. Z,, 



2 Hij 
E 



(12) 



where E$ is the Coulomb barrier energy 56]. Although 
in the standard cross section parameterization, shown 



in equation (|llfl . the second Coulomb penetrability fac- 
tor can be absorbed into the astrophysical S factor be- 
cause, typically, the energy is much smaller than the Q- 
value, here we retain it explicitly. The phase space factor, 
<$> kl (E + Q),is 



$ k i(E + Q) oc J(E + Q)n 3 kl 



(13) 



while the statistical weight factor g k i accounts for the 
multiplicity of the final state 



g u = (2J fe + l)(2J ; + l) 



(14) 



where J k and J\ are the spins of the individual nuclei. 
The reciprocal of the entrance channel velocity is pro- 
portional to <&ij/(nijE). Putting all these together the 
cross section for the reactions i + j <-> k + l is expected 
to be of the form 

(E) = S l3M (E)^-G tJ (E)G kl (E + Q) 

(15) 



x$y(£) <S> M (E + Q). 



where we have introduced Sij^i{E) as an undetermined 
function. This function becomes constant at low energy, 
and is similar to, but not the same as, the astrophysical 
S-factor. 

With the help of equation Ijl5|l we can extract the 
most obvious behavior of any cross section with energy 
to derive Sij^ki(E). Our expectation is that this quantity 
varies slowly though of course the exact details of a reac- 
tion may lead to significant departures. After extracting 
Sij,ki{E) from some known reaction we can then insert 
it into the similar dineutron reaction based on the as- 
sumption that Sij t ki does not change considerably from 
one to the other. The major changes in the cross sec- 
tions will therefore be limited to the considerable effects 
of the Coulomb barrier penetrability and phase space. 
We have examined the validity of this approach by using 
it to predict D(D, p)T from D(D, n) 3 He, and T(D, n) 4 He 
from 3 He(D, p) 4 He. The last two cross section are domi- 
nated by large resonances corresponding to excited states 
of 5 He and 5 Li 59] but the non-resonant pieces have been 
extracted by Chulick et al. |57| allowing us to compare 
the transformation. In both test cases we find the trans- 
formation works reasonably well with an error that is a 
factor of order a few. 

For the 2 n(p, n)D reaction there is no similar deuteron 
reaction with which to compare so instead we used the 
broadly similar 3 He(n, p)T. We have been unable to find 
an analytic expression for the this cross section so we 
interpolated the ENDF-IV evaluated cross section data 
available online 58] and then factored out the expected 
behavior shown in equation (|15|) before replacing it with 



the appropriate terms for 

2 



f p <-> D + n. 

We expect the D + 2 n <-» T + n cross section to be 
similar to D(D,p)T and D(D, n) 3 He up to corrections 
for the Coulomb barrier penetrability and phase space 
factors. Here we have analytic expressions of the S-factor 
to use from Chulick et al. I57I. 
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To estimate the last two reactions, 3 He + 2 n <-> n+ 4 He 
and 3 He + 2 n <-> D + T, one appeals to their similarity 
with T(D, n) 4 He so that one may use the Chulick et al. 
|57| expression and, once again, correct for the change 
in the phase space, Coulomb barrier penetrability etc. 
As mentioned earlier, the T(D, n) 4 He cross section ex- 
hibits a resonance due to an excited state of the 5 He 
nucleus (see [59j for an energy level diagram). The po- 
sition of this same resonance, relative to the 3 He + 2 n 
ground state, depends on the dineutron binding energy 
being subthreshold for £?2 n < 3 MeV. In addition, there 
are further excited states of 5 He that become relevant 
when i?2 n ~ but, as we will show, the effects of the 
dineutron become apparent only when B2 n approaches 
£?d and we have not added them to our cross section. 



D 



This weak reaction is actually the sum of four sub- 



processes: 



D 



"n + v e <-> D 



< E v < A2 nD — m e (16a) 
-A2 nD + m e < E v 



"n + e <-> D + u e , A2 nD + m e < E v 



(16b) 
(16c) 

2 n + e + v e «-> D < E v < - A2 nD - m e (16d) 

where we have denoted by A2 nD the dincutron-deuteron 
mass difference i.e A2 nD = A ap +Br>— Bz n = 3.52 MeV — 
£?2 n . We do not consider those cases where the final states 
are free nucleons. From the limits on the neutrino energy, 
E v , we see that, at most, only three of these reactions can 
be operant at any given value of i?2 n . For these rates we 
use expressions similar to the neutron-proton intercon- 
version rates but with different Q-values and a matrix 
element that takes into account the presence of two neu- 
trons. We also use a pure Gamow- Teller decay between 
the + ground state of the dineutron and the 1 + ground 
state of the deuteron. While there remains some uncer- 
tainty in the rates it has a much smaller impact on final 
abundance yields as compared with the uncertainties in 
the strong interaction rates discussed above. 

There is one interesting quirk that appears when 
A2 nD < m e , which is that an atom of deuterium can cap- 
ture an electron to form a dineutron. Although interest- 
ing this process has not been included in our calculations 
because the amount of atomic deuterium is negligible at 
the temperatures relevant to BBN. We also see that spon- 
taneous deuteron decay can occur when A2 nD < — m e i.e. 
when £?2 n > 4 MeV. Though this possibility remains in- 
teresting we shall not pursue dineutron binding energies 
that permit this reaction to occur. 



no charge. Despite this smallness, it is the only exother- 
mic reaction capable of producing dineutrons when B2 n 
is small. Following Rupak 60j, the lowest order con- 
tribution to this cross section should be something like 

{N T a 2 ® T 2Tj Ny(N T (7 2 <g) T 2 T j {D k + D k )N)E k and 

(N T a 2 <8 T 2 T j N^(N T a 2 (T i ® T 2Tj (D k - D k )N)Bit lkl . 
From examining the n(p, 7)D operators in the cross sec- 
tion from Rupak we estimated that this lowest order con- 
tribution is N LO. We therefore make an order of magni- 
tude estimate for the dineutron cross section by starting 
with the n(p, 7)D (which has a leading order contribu- 
tion) and suppressing it by the appropriate factor. Al- 
though there is considerable error, just as with the two 
3 Hc + 2 n reactions, the effects of the dineutron will only 
become apparent when £?2 n approaches B^> by which time 
this reaction will be of little importance. 

The 2 n(p, 7)T reaction has been included as a failsafe 
mechanism to remove dineutrons since it is exothermic 
for all i?2 n . It is not expected to be an important reaction 
unless we strongly deviate from the nucleoli flow in stan- 
dard BBN because the similar, standard BBN reaction 
D(n, 7)T is also unimportant. We have also estimated 
its cross section from the n(p,7)D reaction in 60] after 
modifying the phase space and spin multiplicity factors 
using the expressions discussed above. 



IV. THE EFFECTS OF A BOUND DINEUTRON 
UPON BBN. 

The chief manner in which the dineutron affects BBN 
is via the Q-values for the reactions and, from the values 
quoted earlier, we can identify three regions of B2 n up to 
the 4 MeV limit we considered. They are: 

• B2 n < B u , 

• B u < B2 n < 3.0 MeV, and 

• 3.0 MeV < B 2n . 

In figure J5J) we plot the fractional change in the pri- 
mordial abundance relative to standard BBN at r\ = 
6.14 x 10~ 10 and the three regions we have identified 
are clearly visible. Up to B-q there is no discernable 
change with £?2 n , over the interval B-q < B2 n < 3.0 MeV 
the deuterium and helium-3 abundances drop and the 
helium-4 abundance rises but after 3.0 MeV the helium- 
3 and helium-4 abundances plateau and the evolution of 
Yd changes noticeably. In what follows we shall explain 
why and how the dineutron influences BBN in each re- 
gion. 



n(n, 7) 2 n and 2 n(p,7)T 

Though n(n, 7) 2 n looks similar to deuteron formation 
via n(p, 7)D this reaction is suppressed because there is 



A. B2 n < -E>d 

For £?2 n <C £?d the dineutron lifetime is of order 1 s so 
while the decay would appear to be a significant drain on 
the dineutron abundance any loss is easily, and rapidly, 
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FIG. 2: The fractional change in the primordial deuterium, 
helium-3 and helium-4 abundances as a function of the dineu- 
tron binding energy. 



deuterons, combining for a total of ~ 1/20. The figure 
also makes it quite clear that as £?2 n approaches Bd the 
dineutron starts to possess a considerable abundance at 
the transition to BBN proper, approximately the loca- 
tion of the deuteron peak. And the figure also shows 
that the position of the peak dineutron abundance is 
unaffected by its binding energy and is coincident with 
the deuteron peak. Whatever the dineutron binding en- 
ergy in this range BBN proper is still initiated by the 
deuteron's departure from NSE and the small amount 
of dineutrons present at that time is rapidly removed. 
When i?2 n <c Bn. their small abundance means that they 
never form a substantial population which could possibly 
influence the predictions of BBN but as B2 n approaches 
Bd their presence at Tbbn becomes important. 



Bo < B 2n < 3.0 MeV 



replaced from the pool of free neutrons and, consequently, 
the dineutron abundance during the second stage of BBN 
follows its NSE value: 

1 Y n Y D { Bi D -B B \ 
= 3 ~~Yp~ 6XP [—f— ) ■ 

Therefore the effect of increasing the dineutron bind- 
ing energy is the same as one expects from NSE, i.e. 
the abundance builds up at higher temperature as the 
binding energy increases. The evolution of the dineutron 
abundance as a function of the temperature for the four 
cases Bi n £ {10 4 eV, 10 6 eV, 2 x 10 6 eV,B D } is shown 
in figure @. In every example in the figure the dineu- 
tron abundance is smaller than the deuteron abundance: 
even when B2 n = B& the NSE abundance is smaller 
than Yd due to both the smaller abundance of free neu- 
trons (Y n /Y p ~ 1/6 — 1/7) and the spin factor, 1/3, of 
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FIG. 3: The evolution of the dineutron abundance as a func- 
tion of the temperature for B2 n £ {10 4 eV, 10 6 eV, 2 x 
10 6 eV,i?D} and the deuteron abundance from standard 
BBN. 



As shown above, the NSE dineutron abundance at 
T ~ 100 keV when £?2 n = Bn. is smaller than Yd by 
roughly a factor of 1/20 so it is not until the dineutron 
binding energy has grown to the point where it is capable 
of reversing the suppression of its NSE abundance by the 
neutron-to-proton ratio and the spin of the deuteron that 
its abundance approaches Yjj at T ~ 100 keV. Roughly 
this occurs when 

exp (*^) ~ 20 (18) 

i.e. when B2 n ~ 2.5 MeV. The evolution of the dineu- 
tron and deuteron abundances during this interval for 
£?2 n is shown in figure @ and confirms that the dineu- 
tron abundance surpasses the deuteron abundance at 
£?2 n ~ 2.5 MeV. The figure also demonstrates a num- 
ber of new features: the position of the dineutron peak 
abundance now moves to higher temperatures as £?2 n in- 
creases, the peak deuterium abundance drops and the 
temperature at which it departs NSE also shifts to higher 
temperatures. These effects ripple through to the triton 
and helion as shown in figure JSJ). As B2 n increases tri- 
tium and helium-3 depart NSE at lower temperatures: 
the change isn't dramatic, by B2 n — 3 MeV the two 
A = 3 nuclei depart at T ~ 170 keV rather than at 
around 200 keV in standard BBN but the sensitivity 
of the NSE abundance to the temperature means that 
the abundance of hclium-3 and tritium after this point is 
roughly 2 — 3 orders of magnitude larger. The figure also 
shows that the movement of the tritium peak parallels 
that for dineutrons by shifting to higher temperatures as 
B2 n increases. 

The movements in the abundances of the intermediary 
nuclei are all attributable to changes in the flow of the nu- 
clei through the reaction network. For £?2 n > Bd, the re- 
action D(n, p) 2 n is now exothermic and a detailed exami- 
nation of the nuclear flow at B2 n = 2.6 MeV confirms this 
to be the source of dineutrons. The flow also indicates 
the dineutrons are then chiefly converted to tritons via 
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FIG. 4: The evolution of the deuterium (top panel) and dineu- 
tron (bottom panel) abundances as a function of the tempera- 
ture for i?2 n £ {Bo (dotted line), 2.6 MeV (dot-dashed line), 
3MeV(solid line)}. 



FIG. 5: The evolution of the tritium (top panel) and helium- 
3 (bottom panel) abundances as a function of the tempera- 
ture for £?2 n G {Bo (dotted line), 2.6 MeV (dot-dashed line), 
3 MeV (solid line) }. 



the reaction 2 n(D, n)T. Both reactions lead to a disrup- 
tion of the usual mechanisms that lead to BBN proper. 
In standard BBN the transition to BBN proper was due 
to the removal of deuterons via the two D + D processes 
but now, because beyond B~2 n ~ 2.5 MeV the NSE abun- 
dance of dineutrons is larger than that of deuterons and 
because free neutrons are so plentiful during the second 
stage of BBN, the departure occurs earlier. Remarkably 
both D(D, n) 3 He and D(D,p)T are now of little conse- 
quence. The shift in the deuteron's NSE departure tem- 
perature is not large, figure @J demonstrates this, but, 
as with T and 3 He, the deuteron's NSE abundance is 
very sensitive to the temperature. The 2 n(D, n)T reac- 
tion can also provide sufficient tritium to keep its abun- 
dance in equilibrium until a slightly lower temperature. 
The switch in the mechanism that initiates BBN proper 
explains many of the features seen in figures (@J and JSJ. 

From the nuclear flow we also find that helium-3 no 
longer plays an important role in the formation of helium- 
4. Some helium-3 is still formed, via D(D, n) 3 He, but this 
source is suppressed due to the lack of deuterons; what 
little helium-3 is created is processed to tritium by the 
usual 3 He(n, p)T though 2 He( 2 n, n) 4 He does play a role. 



The reduced significance of helium-3 is not reflected in 
figure ®: one must remember that the net rate of for- 
mation for the intermediary D, T and 3 He is the small 
difference between production and destruction and does 
not necessarily indicate the true amount of nucleons flow- 
ing through them. Even in standard BBN the evolution 
of helium-3 does not resemble that of tritium because 
3 He(n,p)T is so rapid and one can only see the signifi- 
cance of the helium-3 nucleus by examining the individ- 
ual reaction rates |5l) . 

Finally, the formation of helium-4 still occurs via 
T(D, n) 4 He and, due to the earlier initiation of BBN 
proper, its final (primordial) abundance is enhanced. The 
reaction network is modified and the bulk of the nucle- 
ons now flow through a network resembling that shown 
in figure I©- 



By B2 n = 3 MeV figure @ shows that the peak abun- 
dance of deuterons has dropped by two orders of magni- 
tude from that in standard BBN and BBN proper begins 
at an even higher temperature. The nuclear flow now in- 
dicates that tritium formation is increasingly dominated 
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FIG. 6: The schematic flow of nuclei through the reaction 
network from deuterium to tritium at £>2 n = 2.6 MeV. The 
flow does not include helium-3. A little of this nucleus is 
still produced via D(D, n) 3 He but the suppressed deuteron 
abundance means that the amount is substantially reduced 
compared to standard BBN. 



by 2 n(p, 7)T rather than 2 n(D, n)T and, interestingly, the 
dominant source of helium-3 is the mildly endothermic 
T(p, n) 3 He at higher temperatures with a changeover to 
T(D, 2 n) 3 He as the Universe cools. This flow, T(p, n) 3 He, 
is opposite to that in standard BBN, 3 He(n, p)T. 

This second region of Bi n is where the presence of 
dineutrons really becomes manifest. At i?2 n = Bd the 
dineutron is just starting to influence BBN, by £?2 n = 
3 MeV it has considerably altered the reaction network 
that take the free neutrons to helium-4 and led to a signif- 
icant change in the mechanism that leads to BBN proper. 



3.0 MeV < Bi 



only through the exit channel phase space which varies as 
-\Z_Bt — B2 n . With an abundance controlled by D and 
a destruction mechanism that varies only weakly with 
B-2 n the height and temperature of the peak dineutron 
abundance is, essentially, static. Deuterium departure 
from NSE occurs when its abundance opens the D(n, p) 2 n 
drain and, again, B2 n only enters weakly through the exit 
channel phase space. 

A new feature emerges in this third region of B2 n . At 
B2 n > 3.0 MeV, the Q-value for the reaction 3 He( 2 n, D)T 
is now negative and the importance of the helion as 
an intermediary in the formation of helium-4 rebounds, 
though not to the level in standard BBN. The source of 
helions is principally the mildly endothermic T(p, n) 3 He 
but switches to T(D, 2 n) 3 He as the Universe cools. 

Finally, helium-4 is now chiefly formed by both 
T(D, n) 4 He and T(T, 2n) 4 He in almost equal propor- 
tions with minor contributions from T(p,7) 4 He and 
3 He( 2 n, n) 4 He. The network has changed slightly and 
the new flow at B2 n = 3.5 MeV is illustrated in figure 
0. 
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FIG. 7: The diagrammatic flow of nuclei through the reaction 
network from tritium onwards at E>2„ = 3.5 MeV. 



The excitement seen when Bn < £>2 n < 3 MeV 
has largely played out as we enter the third domain of 
3 MeV < £?2 n and the shifts in the evolution of the abun- 
dances slows. In this range, between 3 MeV and 4 MeV, 
there are no, new, major changes in the nuclear flow be- 
cause the Q-values of the most important reactions are 
now all ~ MeV. Nevertheless, a detailed study of the 
reaction network does show small differences and we dis- 
cuss those here. 

Although further increases in the dineutron binding 
energy are not reflected in the position and amplitude of 
the peak dineutron abundance, the temperature at which 
the dineutron departs NSE shifts to higher values with 
£?2 n . In contrast, the reduction in the amplitude of the 
peak deuteron abundance seen in figure @ becomes less 
dramatic and its NSE departure temperature has all but 
ceased to move as B2 n increases. 

The behavior of the evolution of these two A = 2 nu- 
clei reflect the fact that dineutrons are primarily cre- 
ated from deuterons and that the dominant dineutron 
destruction mechanism has switched from 2 n(D, n)T to 
2 n(p,7)T. Dineutron NSE departure occurs because of 
insufficient production from the small D abundance. The 
2 n abundance at lower temperatures is thus controlled by 
the deuteron. The reaction 2 n(p, 7)T is affected by B2 n 



V. ARE THESE PREDICTIONS ROBUST? 

While we have strived to estimate the various dineu- 
tron cross sections by basing them upon simple physical 
assumptions, it is a worthwhile task to see how robust 
the effects we have described are. If we know (or assume) 
the distribution for the error, either in the cross-section 
or the thermally averaged rate, then one can compute the 
error matrix as in Fiorentini et al. |62| and Cuoco et al. 
|6l| though this approach cannot recover higher moments 
of the abundance distribution. An alternative is to con- 
struct the distribution of the abundances by sampling 
the distribution of the errors in a Monte Carlo analy- 
sis. This technique was used by Krauss & Romanelli 
1(53], Smith, Kawano & Malaney [Hi]], Krauss & Kernan 
|64j and, more recently, Nollett and Buries [6^ and is 
the technique we shall use. To this end we introduced 
random multiplicative factors for all our dineutron reac- 
tion cross sections with the exception of 2 n <-* D, the 
most reliably estimated. These factors were chosen from 
a probability distribution limited to the range between 
1/5 and 5 and weighted such that there was equal prob- 
ability either side of 1. The baryon-to-photon ratio was 
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fixed at r\ = 6.14 x 10 . In figure 1(8} we plot the root 
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FIG. 8: The root mean square deviation divided by the mean 
plus the ±3<J error for deuterium and helium-4 as a function 
of the dineutron binding energy. Below E>2 n ~ 2 MeV nei- 
ther nucleus exhibits any spread due to the uncertainty in 
the dineutron reaction rates. Above this value the spread 
in the deuterium results becomes large while the results for 
helium remain smaller than ~ 2% for all value of Bi n . 



mean square deviations for deuterium and helium-4 and 
the three regimes for the dineutron binding energy are 
clearly visible in both curves. The increasing spread for 
both reflects the increasing dominance of dineutron reac- 
tions in the formation of helium-4. The figure shows that 
up to Bi n r~i 2 MeV the errors in the reactions involving 
dineutrons do not introduce any error into the predicted 
abundances, illustrating again the insignificance of dineu- 
trons in BBN when their binding energy is smaller than 
this value. Above £?2 n ~ 2 MeV the curves begin to de- 
viate from zero but note two important features: first, 
the spread in helium-4 is small for all values of Bz n , and 
second, the spread in the predicted abundance of deu- 
terium does not exceed ±100%. The small spread in the 
helium-4 curves indicate that we can reliably predict the 
abundance of hclium-4 even if dineutrons have signifi- 
cantly impacted BBN, while the fact that the deuterium 
curves do not exceed 100% shows that at least the direc- 
tion of the change is known even if the exact abundance 
is not. We have not shown the rms spread for helium-3 
because it exceeded ±100%. 

Note also that the range in the abundances is sig- 
nificantly smaller than the range we permitted for the 
random factors. This may seem remarkable since the 
primordial abundance of an intermediary nuclei, such 
as deuterium, involves a fierce competition between its 
production and destruction and hence the adopted cross 
sections. While dineutrons may have led to significant 
departures in the flow of nuclei compared to standard 
BBN, large variations in the dineutron reaction rates do 
not translate into equally large spreads in the abundances 
of the intermediaries. 

We also found that there was a significant anti- 
correlations in the results as shown in figure (JSJ. For 
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FIG. 9: The correlation between the deuterium abundance 
and helium-4 mass fraction, Pdy, and its ±3cr error, as a func- 
tion of the dineutron binding energy. Below B2 n ~ 2 MeV 
the correlation is close to —100% while after £?2 n ~ 3 MeV 
the correlation has softened to ~ —80%. The large peak at 
Bi a ~ 2 MeV corresponds to the point where the change in 
the flow of nuclei through the reaction network occurs. 



£?2 n < 2 MeV the correlation coefficient was almost ex- 
actly -1 while over the interval 2.5 MeV < B2 n < 3 MeV 
the anti-correlation softened slightly to ~ —0.8. As a con- 
sequence of this correlation the covariance matrix, Vr, 
describing the error in the predictions arising from the 
uncertainties in the dineutron's reaction rates, contains 
off-diagonal pieces. 
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FIG. 10: The fractional difference between the mean from 
the Monte Carlo simulation and the result with no random 
factors. The solid line is for helium-4, the dotted is for deu- 
terium. 



We have delayed until last the difference we find in the 
mean from the sample and the primordial abundances we 
derive with no random factors as shown in figure (J2J . The 
fact that a discrepancy arises is not surprising given the 
non-linearity of BBN. The errors (which are functions of 
£?2 n ) for helium-4 is < 0.5% but for deuterium it reaches 
< 20% and helium-3 fared even worse with a discrepancy 
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between the mean and no random factors approaching 
30%. The fractional difference between the mean and 
the result with no random factors are shown in figure 
(|10|l for deuterium and helium-4. Although seemingly 
large, this error is smaller than the statistical fluctuations 
for all three nuclei (seen in figure JBJ for deuterium and 
helium-4) but not significantly so and we must include it 
in the total error for the predictions as a systematic. The 
covariance matrix for the systematic error is denoted by 
V s . 

With means, Y, and variance, Vr = Vr + Vs, that are 
functions of Bi n we approximate the distributions in the 
predictions as Gaussians 



P(Y\B2 n ) = 



x exp 



(Y - Y) T VjT 1 (Y - Y) 



(19) 



We have checked the validity of this approximation to the 
spread in the results by performing a Kolmogorov test at 
each value of Bi n we set. The spread in the predictions 
for helium-3, which we have not shown, did not pass this 
due to a significant kurtosis. 



VI. A DEGENERACY WITH r?? 

So far we have restricted our discussions to a fixed 
value of the baryon-to-photon ratio, r\, but if we allow 
this parameter to also vary we may well end up with a 
degeneracy that makes it difficult to establish the pres- 
ence of a bound dineutron. To determine if this occurs, 
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FIG. 11: Helium-4 iso-mass fraction (solid) and Deuterium 
iso-abundance (dashed) contours in the i?2 n — n plane. From 
top to bottom, the helium-4 contours are Y = 0.248 and Y = 
0.244 and the deuterium abundances are T)/H — 1.8 x 10~ 5 , 
D/H = 2.6 x 10~ 5 and D/H = 3.4 x 10~ 5 . 

we show in figure (|11|) iso-abundance and iso-mass frac- 
tion contours for deuterium and helium respectively as a 
function of Bi a and r\. If we examine the most robust pre- 
diction of our modified BBN, the increase in Y , then we 



can quite easily mask this effect by lowering r\ as shown 
in the figure. The decrease in r\ required to offset the 
increase in the primordial mass fraction as Bi n increases 
is considerable so that above Bi Yl ~ 2.5 MeV all values 
of rj in the 3 x 10~ 10 to 8 x 10~ 10 range plotted yield 
a helium-4 mass fraction above 0.248. For deuterium, 
lower values of rj lead to increases in the final abundance 
so that they too can mitigate the decrease in Yr> as £?2 n 
increases. But the figure shows that the decrease in rj 
required for deuterium is nowhere near as large as that 
required for hclium-4, so, while there is an anticorrelation 
of i?2 n and rj for both helium-4 and deuterium, there is 
no degeneracy for both simultaneously. If we considered 
each species separately then the degeneracy could, in- 
stead, be broken by using the CMB since lowering tjbbn 
will lead to a discrepancy with ijcm b ■ 



VII. AN UPPER LIMIT FOR B 2n 

It's finally time to derive an upper limit to £?2 n based 
on the compatibility with observations. The primordial 
abundance, D/H, of deuterium is taken to be D/H = 
(2.6 ± 0.4) x 10~ 5 [3] while we will consider both the 
Olive, Steigman and Walker 68] value of Yosw — 0.238± 
0.005 and the Izotov and Thuan value of Y IT = 0.244 ± 
0.002 [69ll70| for the helium-4 mass fraction Y. The exact 
primordial abundances remain a topic of debate with two, 
largely incompatible, determinations for the helium mass 
fraction |6^, [7(3, EH IIS an d excessive scatter in the 
measurements of deuterium [lil 0] but these two nuclei 
still represent the best probes of BBN. 

The error in the observed hclium-4 mass fractions are 
2% for Yosw, 1% for Yjt which compares well with the 
spread in the predicted mass fraction plotted in figure 
©■ We can integrate over the possible values for the 
prediction, at a fixed r\ and i?2 n , using the distributions 
found earlier, and find 



C(r,,B2 n \Y) 



1 



x exp 



(Y - Y) T V -1 (Y - Y) 



(20) 



where Y denotes the vector whose elements are the ob- 
servations and Y the vector whose elements are the pre- 
dictions and V is the covariance matrix - the sum of the 
(diagonal) covariance matrix for the observations, Vo, 
and the two covariance matrices Vr and Vs- Contours 
of the likelihood are shown in figure (|12fl which shows 
that the use of deuterium and helium-4 breaks the de- 
generacy seen in each separately, and, indeed, the limits 
to the dineutron binding energy are independent of the 
baryon-to-photon ratio. 

Marginalizing over the baryon-to-photon ratio, r\, we 
obtain the results shown in figure \Y&\. The two upper 
limits for £?2 n are not greatly influenced by the different 
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FIG. 12: The 95% and 99% confidence contours using the 
Olive, Steigman and Walker (solid curves) and Izotov and 
Thuan (dashed curves) helium-4 mass fractions observations 
and the Barger et al. [l4| primordial deuterium abundance. 
The best fit points, both at £>2 n = 0, are denoted by the 
triangle for OSW, the square for IT. 



primordial helium-4 mass fractions found in the litera- 
ture. 
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FIG. 13: The marginalized likelihood function for B2 U using 
the Olive, Steigman and Walker (solid curve) and Izotov and 
Thuan (dashed curve) helium-4 mass fractions observations 
and the Barger et al. 0] primordial deuterium abundance. 



VIII. SUMMARY AND CONCLUSIONS 

We have examined the role that a stable dineutron may 
play in BBN as a function of its binding energy, Bi n , up 
to 4 MeV. We have estimated the important new re- 
actions that enter into the BBN reaction network and 
examined, in detail, the change in the nucleon flow. We 
find that the range < i?2 n < 4 MeV can be subdi- 
vided into three regions: £?2 n < £>d, Bp < B2 n < 3 MeV 
and 3 MeV < B2 n . The boundaries are due to the sign 



inversion of the 2 n(p, n)D and 3 He( 2 n, D)T Q-values re- 
spectively. Below By> the dineutron has little effect upon 
BBN but as we increased Bi n beyond this value it began 
to interfere with the nucleon flow between deuterons and 
tritium. The hclium-4 mass fraction increased by ~ 10% 
and the deuterium abundance dropped by ~ 40%. Above 
3 MeV the nucleon flow settled into a new pattern, the 
helium-4 mass fraction plateaued to new level and the 
rate at which the deuterium abundance decreases with 
£?2 n slowed. A small change to the nucleon flow from 
tritium to helium-4 was seen. 

We estimated the error in the predictions by sampling 
the distribution of abundances when the dineutron reac- 
tions were multiplied by random factors and found that 
the prediction of an increase in the helium-4 mass fraction 
and decrease in deuteron abundance were reliable. The 
degeneracy between Bi^ and the baryon-to-photon ratio, 
■q, that occurs for Y and D/ 1 H separately was broken 
when both were considered simultaneously. We then con- 
structed the 2-D likelihood function £(r?, Bi n \Y, D I H) by 
using both the OSW [H and IT [UlSi helium-4 mass 
fractions and the Barger et al. [l4T | deuterium abundance 
and then marginalized over r\ to derive the likelihood dis- 
tribution for B2 n . We found that dineutron binding en- 
ergies above ~ 2.5 MeV could not be accommodated by 
BBN within both allowed ranges of Y 

Throughout our calculation we have only permitted 
the variation of the dineutron binding energy. In real- 
ity, whatever the source of the stability of the dineutron, 
naively the binding energies of the other nuclei should 
also change. The range in B-2 n we investigated is much 
larger than the range in £>d that Kncller & McLaugh- 
lin permitted but we notice that the increase in Y and 
decrease in D/H seen there when Bd increased is also 
mimicked by the dineutron. Thus if Bd and i?2 n increase 
in tandem then the stability of the dineutron cannot re- 
verse the increase in Y whatever the relative magnitudes 
of the two binding energies and even though dincutrons 
may alter the nucleon flow. If these two binding ener- 
gies change in opposite senses then the situation is more 
confused: the effect of a decrease in Bd is the immediate 
decrease in the helium-4 mass fraction which may not 
be reversed by the presence of a stable dineutron unless 
£?2 n > Bd- A better understanding of the net effect, 
and a more concrete limit to the permitted variation of 
both the deuteron and dineutron binding energies, could 
be derived if the relationship between them were known 
and we were not forced to have to treat each as inde- 
pendent. Furthermore, if such a calculation shows that 
the predicted helium-4 mass fraction can be consistent 
with observations, without changing the compatibility of 
Vbbn and t]cmb and the deuterium abundance greatly, 
then an examination of the effects upon the primordial 
yield of lithium-7, which, in standard BBN, is overpro- 
duced theoretically, would be in order. Although we have 
not included any dineutron reaction involving nuclei with 
mass above A = 4, one can speculate that the omitted re- 
action 7 Be( 2 n, ncv) 4 He could play a significant role since 
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beryllium-7 (before it decays to 7 Li) is the chief compo- 
nent of the primordial A = 7 isobar yield at 77 ~ 6x 10 -10 . 

In this paper we have shown that the dineutron can be- 
come bound at a level up to that of the deuteron without 
disrupting the standard nuclear flow in BBN or signifi- 
cantly altering predicting BBN abundance yields. Be- 
yond that, changes to the nuclear flow and to predicted 
abundance yields appear. It is important to note that we 
have allowed only one parameter, the dineutron binding 
energy, to vary. Further work on the interdependence 
of cross sections and binding energies in nuclear theory 
would be required to reduce the errors presented here and 



to make a more concrete connection with the underlying 
fundamental constants. 
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